Quantum entanglement and Bell inequalities in Heisenberg spin chains 
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We show that in one-dimensional isotropic Heisenberg model two-qubit thermal entanglement and 
maximal violation of Bell inequalities are directly related with a thermodynamical state function, i.e., 
the internal energy. Therefore they are completely determined by the partition function, the central 
object of thermodynamics. For ferromagnetic ring we prove that there is no thermal entanglement 
at any temperature. Explicit relations between the concurrence and the measure of maximal Bell 
inequality violation are given. 



PACS numbers: 03.65.Ud, 75.10.Jm 

Over the past few years much effort has been put 
into studying the entanglement of quantum systems both 
qualitatively and quantitatively. Entangled states consti- 
tute indeed a valuable resource in quantum information 
processing jjj. Entanglement in systems of interacting 
spins p Mfll a s well as in systems of indistinguishable par- 
ticles |[7Hll| has been investigated. In particular entan- 
glement in both the ground state |^|| and thermal state 
of a spin- 1/2 Heisenberg spin chain have been ana- 
lyzed in the literature. The intriguing issue of the relation 
between entanglement and quantum phase transition Jl2] 
have been addressed in a few quite recent papers [ ]T3| . |14[ . 

From a conceptual perspective one can say that these 
investigations aim to provide a bridge between quan- 
tum information theory and theoretical condensed mat- 
ter physics. This is done by considering thermal equi- 
librium in a canonical ensemble. In this situation the 
system state is given by the Gibb's density operator 
p T = cxp (—H/kT) /Z, where Z =tr[exp (-H/kT)\ is the 
partition function, H the system Hamiltonian, k is Boltz- 
mann's constant which we henceforth will take equal to 
1, and T the temperature. As p(T) represents a ther- 
mal state, the entanglement in the state is called ther- 
mal entanglement p|. It is important to stress that, al- 
though the central object of statistical physics, the par- 
tition function, is determined by the eigenvalues of H 
only, thermal entanglement properties generally require 
in addition the knowledge of the energy eigenstates. 

On the other hand the violation of Bell inequalities ]lq] 
was considered as a means of determining whether there 
is entanglement. In 1989 Werner |l(| demonstrated that 
that there exists states which are entangled but do not 
violate any Bell-type inequality, i.e., not all entangled 



states violate a Bell inequality. Further studies [ |L7| 18| 
showed that the maximal violation of a Bell inequality 
does not behave monotonously under local operation and 
classical communications. So Bell violations can only be 
considered an entanglement witness Jl9| . More recently 
the relation between Bell inequalities and the usefulness 
for quantum key distribution and quantum secret sharing 
have been clarified pfj[ . Moreover it has been showed 
that even multipartite bound entangled states can violate 
Bell inequalities |21[. 



In this paper we shall study the Heisenberg qubit 
chains by exploring further connections between entan- 
glement and other relevant physical quantities. We will 
study the relations between the concurrence, an entan- 
glement measure for two qubits, with thermodynamic po- 
tentials such as the internal energy U and magnetization 
M. The latter quantities are defined as 
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where (3 = 1/T and B is a magnetic field. We shall 
establish a remarkably simple relation between the con- 
currence for the thermal state and the internal energy. A 
direct relation between the maximal violation of the Bell 
inequality and the internal energy will be given as well. 

We first briefly review the definition of concurrence 
p2| . Let pi2 be the density matrix of a pair of qubits 1 
and 2. The density matrix can be either pure or mixed. 
The concurrence corresponding to the density matrix is 
defined as 



C = max {Ai — A2 — A3 — A4, 0} . 



(2) 



where the quantities Ai > A2 > A3 > A4 are the square 
roots of the eigenvalues of the operator 



(3) 



The nonzero concurrence implies that the qubits 1 and 
2 are entangled. The concurrence C12 = corresponds 
to an uncntangled state and C12 = 1 corresponds to a 
maximally entangled state. 

The Heisenberg model and pairwise entanglement. We 
consider a physical model of a chain of N qubits, namely, 
a chain of spin-i particles in which neighboring particles 
interact via the anisotropic Heisenberg Hamiltonian H = 
H(A, B), with a magnetic field 
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N 



H 
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(4) 
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where cr, = (a ix , o~ iy , o~ iz ) is the vector of Pauli matrices 
and J is the exchange constants. The positive and neg- 
ative J correspond to the antiferromagnetic (AFM) and 
ferromagnetic (FM) case, respectively. 
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In order to enact translational invariance for finite 
N we assume periodic boundary conditions, i.e., N + 
1 = 1, turning the chain into a ring. Therefore H 
becomes invariant under cyclic shifts generated by the 
right shift operator T. The latter being defined by its 
action on the product basis, T \mx, TOjv— i> ^at) = 
\niff, mi, rniv-i)) where m, = (1) represents the 
state of spin up (down). Another fact is that that 
[H, S z ] = 0, which guarantees that reduced density ma- 
trix pi2 of two nearest-neighbor qubits, say qubit 1 and 
2, for the thermal state has the form [SI 



P12 



Wi z 
Z W 2 



(5) 



in the standard basis {|00), |01), 1 10) , |11». Here S a = 
J2iLi °~a/2 (a = x, y, z) are the collective spin operators. 
The reduced density matrix is directly related to various 
correlation functions Gap = (ai a a 2 p) =tx(ai a (T20 pr) ■ 
Then matrix elements can be written in terms of the cor- 
relation functions and the magnetization per site M — 
M/N as 



•± = i(l±2M + G zz ), 



(6) 



In the deriving of above equation, we have used the trans- 
lation invariance of the Hamiltonian. From the above 
equation we find the following relations which will be use- 



(l + G„),He(z) 



ful later u + — u = M, u + + u - 
\{G XX + G yy ). All the information needed is contained in 
the reduced density matrix, from which the concurrence 
quantifying the entanglement is readily obtained as 



C = 2max[0, \z\ - Vi 



+■,,- 



)}, 



(J) 



which can be expressed in terms of the correlation func- 
tions and the magnetization as seen from Eq.(^J). 

Isotropic Heisenberg model. Now we consider the 
Heisenberg XXX model described by the Hamiltonian 
H(1,0). We first notice that H admits a continuous 
SU(2) group of symmetries. This can be easily checked 
in that H commutes with the Lie-algebra generators S a , 
In particular the isotropic Hamiltonian also commutes 
with the operators 



Q a = <J® N (a = x,y,z) 



(8) 



that generate Z2 sub-groups of SU(2). The rich rota- 
tional symmetry structure of the XXX model along with 
translational invariance will play a key role in our study 
of the ring entanglement. 

Form the Z 2 symmetry immediately follows 
Proposition 1. For any temperature the magnetization 
M is vanishing and the nondiagonal clement z is real. 



Proof: By definition one has M 



NZ-hi(a lz e-P H ) = NZ-Hi{o-® N o- u o-® N e-P H ) 



-NZ-Ho-u 



-M. Moreover z 



{0-1+0-2-/ 



((72+C1-) = Z*. □ 

In passing we notice that the reality of the correlation 
function (a^crj) can be proven in a more general context 
possibly with ^-symmetry broken. Indeed it is easy to 
check that it is sufficient that there exist a symmetry S of 
the Hamiltonian swapping the j-th with the i-th lattice 
site i.e., S S = af. This is the case, for example, 
for the XXX on a ring with S given by a geometrical 
symmetry of the lattice, e.g., a reflection. 

In view of Prop. 1 the concurrence has the follow- 
ing simplified form, C = 2max[0, \z\ — u + ]. Furthermore, 
from the relations for u, v, and z established above we 
find that the concurrence is completely determined by 
the correlation functions G aa {a = x,y, z), i.e., 



C = imax[0,|G a 



Gyy\ G z 



-!)]• 



(9) 



This result is due to the Z2 symmetry generated by 
crf N . By using the global SU(2) symmetry and the 
translation invariance it is starightforward to check that 
G xx = Gyy = G zz and G xx = U/(3JN). Then we get 

Proposition 2 For the isotropic Heiserberg model the 
concurrence of nearest-neighbor qubits is directly ob- 
tained from the internal energy, i.e., 



C 



max[0,2|J7/J| - U/J— 3)]. 



(10) 



In Ref. H it has been shown that in the Heiserberg 
even-site antiferromagnetic (J = 1) XXX model g] 
there exists a direct relation between the pairwise en- 
tanglement and the ground state (GS) energy Eos'- 
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C = - max (0, ~E GS /N - 1) 



(11) 



The key step for the proof of this relation is to prove 
that the nondiagonal element z < for the ground state. 
Note that this relation is obtained for the case of anti- 
ferromagnetic and even-number qubits. One can guess 
that the nondiagonal element z is also negative for the 
thermal state of the antiferromagnetic model. In fact we 
can prove a more general result given by 

Proposition 3 The internal energy is always negative 
and the off-diagonal element z is negative (positive) for 
the antiferromagnetic (ferromagnetic) rings. 

Proof: First, from the tracelessness of the Pauli opera- 
tors, it is immediate to check that in the limit of T — * 00 
one has {7 = 0. The internal energy is known to be a 
non-decreasing function of the temperature (dU/dT = 
(AH) 2 /T 2 > 0) U is negative. Since we know that 
z = U/(6JN), we conclude that z < for the antifer- 
romagnetic case and z > for the ferromagnetic case 
□ 

The combination of the Prop. 2 and 3 gives rise to 
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Theorem 1 The concurrence of the nearest-neighbor 
qubits in the Heisenberg model is given by 



a max t ' Sn ~ l \ 



for AFM, 
for FM. 



(12) 



In the limit T — *■ in the antiferromagnetic even- 
site model, the thermal state become the nondegenerate 
ground state. As we expected Eq.(|l^) reduces to Eq.([ll|) 
in this limit. Theorem 1 establishes a direct relation be- 
tween the concurrence and a macroscopic thermodynam- 
ical function, the internal energy. Then the entangle- 
ment is uniquely determined by the partition function of 
the system. Notice that Theorem 1 has no restriction to 
J > and N even (like in Ref §]). Eqs. @ clearly 
show that in AFM (FM) rings the concurrence is a not 
increasing (decreasing) function of the temperature. 

From the obvious fact 

\{a z ®a z )\ = Z- l \^e-P Ei {i\a z ®a z \i)\ < (13) 



-0E t 



one has that — 1 < U/3JN = G zz < 1. Therefore we 
obtain a general result which is applicable for both even 
and odd number of qubits. 

Corollary 1 At any temperature there is no thermal en- 
tanglement between two qubits in the FM XXX Heisen- 
berg rings. 

Now we discuss the case of AFM. When the temper- 
ature increases the internal energy will increase but the 
concurrence decrease. When the internal energy arrives 
at —N J the concurrence becomes zero. 

The temperature T c at which the concurrence vanishes 
is called the threshold temperature. In a short summary 
we have 

Corollary 2 At temperatures lower than T c there exists 
thermal entanglement between two qubits in the AFM 
XXX Heisenberg rings. 

T c is determined by the equation u{N) := 
U(T c )/(—NJ) — 1. From the numerical evidences of 
Refs. Q and |^3| we conjecture that u(N) is a non- 
increasing (non-decreasing) function of N for N even 
(odd). Therefore since it has been estimated u(oo) > 1 || 
and u(3) < l,u(5) > 1 |23| one obtains that the ground 
state of the AFM Heisenberg is always entangled except 
for the case N = 3. 

Now we give the examples of 2 and 3 qubit for the il- 
lustration of our general results. For two-qubit model the 
partition function is given by Z = 3e _2 ^ J + e 6 ^ J , then 
the internal energy follows 6J(e~ 2/3J — e 6 ^ J )/Z. From 
Eq.(|l2"|) the concurrence is found to be zero for the FM 
rings and C = max{0, [exp(8/3J) - 3]/ [exp(8^J) + 3]} 
for the AFM Q. The threshold temperature is then 
determined by exp(8/5J) = 3 and it is given by T c = 
8J/ln3. For three-qubit model the partition function is 
given by 8 cosh(3/3J), and then the internal energy is 



— 3Jtanh(3/3J), from which and Eq.([12|) we can find that 
the concurrence is zero for the FM rings as we expected 
and C = max{0, tanh(3/3J) — l}/2. Hence we recover the 
result that there is no pairwise thermal entanglement for 
the three-qubit Heisenberg ring || . Next we discuss Bell 
inequality. 

Bell inequality. The most commonly discussed Bell- 
inequality is the CHSH inequality |J]|§. The CHSH 
operator (a, a', b, b' are unit vectors) reads 

B = a-a® (b + 9) ■ a + a' ■ a <g> (b - 9) ■ a. (14) 
In the above notation, the Bell inequality reads 



(B) = tr(pB) < 2, 



(15) 



where p is an arbitrary two-qubit state. The maximal 
amount of Bell violation of a state p is given by |2a] 



B = 2Vu + u, 



(16) 



where u and u are the two largest eigenvalues of T p T^. 
The matrix T p is determined completely by the correla- 
tion functions being a 3 x 3 matrix whose elements are 
(T p ) nm =tr(p 12 o- n ® a m ). Here a x = a x .a 2 = <J y , and 
(73 = a z are the usual Pauli matrices. We call the quan- 
tity B the maximal violation measure, which indicates 
the Bell violation when B > 2 and the maximal viola- 
tion when B = 2y/2. The violation measure is a function 
of the correlation functions and does not depend on the 
magnetization. 

For the isotropic Heisenberg model the matrix T p is 
easily obtained as di&g[G xx , G xx , G xx ] and the violation 
measure becomes 



B = 2V2\G X 



2^ For AFM, 
For FM. 



3JN 
2V2U 
3JN 



(17) 



As we expected, the violation measure is completely de- 
termined by the internal energy just as the concurrence. 

From Eqs. (g) and @ we arrive at an explicit relation 
between the concurrence and the violation measure 



C 



imax[0,^ 



1] 
1] 



For AFM, 
For FM. 



(18) 



From the above relation we know that the concurrence 
is larger than zero when the Bell inequality is violated. 
When 2%/2/3 < B < 2 for the AFM rings the state is 
entangled, but the Bell inequality is not violated. The 
result is general for arbitrary N. The maximal value of B 
is 2v2, therefore the concurrence is zero for the FM rings. 
Next we consider the general model with anisotropy and 
magnetic fields. 

General models including anisotropy and magnetic 
fields. For A ^ 1 one obtains an anisotropic model which 
has no longer a global SU(2)- symmetry. The correspond- 
ing Hamiltonian i?(A, 0) still commutes with the z com- 
ponent of the total spin. So the concurrence is given by 
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Eq.(^|). For this model we have the following relations 

U I J = G XX + Gyy " 

nal energy per site. 
Eq.(§) gives 



AG ZZ , where U = U /N is the inter- 
Thc combination of the relation and 



C 



- max[0, \U/J— AG Z 



G z 



1)]. 



(19) 



Then if we calculate the concurrence, we need to know 
the correlation function G zz and the partition function. 
The partition function itself is not sufficient for determin- 
ing the entanglement. In particular for XY model, i.e., 
A = 0, the concurrence reduces to C — \ max[0, \U/J\ — 
G zz — 1)1 . We still need to know the correlation function 
G zz [£6| to calculate the concurrence. 

Now we consider the XXX model with a magnetic 
field, i.e., the Hamiltonian H(l, B). Now the magnetiza- 
tion is no long zero and hence u + ^ u~ . For this model 



we have the relation U / J 

- u 
G 



Gyy + G Z 



BM 



Gxx 

4 Re(z) + 2(u + + u~) - 1 + BM. Here we have used the 
relation G xx + G yy + G zz = 4Rc(z) + 2(u + + u~) - 1. 
Since u + — u~ = M, we can express u and v in terms 
of z ,U I J, and M. Finally the concurrence is C = 
2 max{0, \z\ - \[{U / J - 4Re(z) + 1 - BM) 2 - AM 2 ] 1 / 2 }. 
To calculate the entanglement we need to know, beside 
the partition function, the non-diagonal element z. 

Comments. Our discussions above are applicable to 
more general Heisenberg Hamiltonians such as 



H 



(20) 



where J"j((x = x, y, z) are arbitrary exchange constants. 
For this model we still have the Z 2 symmetry [H, a® N ] = 
0. Therefore the concurrence for the two qubits i and 
j of the thermal state is given, with an extra i,j de- 
pendence, by Eq.(^|). If Jfj = JV. = J[- } , then it 
is easy to check that the Hamiltonian H has a global 
S'J7(2)-symmetry. Therefore the concurrence becomes 

1)], which is determined 



C = imax[0,2|G, z | - G z 
solely by the correlation function G zz . 

Conclusions. In this paper we have discussed ther- 
mal entanglement and Bell inequality violation in the 
Heisenberg qubit rings. In the isotropic case we found 
that there exists a simple relation between the pairwise 
entanglement for the thermal Gibb's state and the in- 
ternal energy. This result is noteworthy and somewhat 
surprising in that it allows to directly connect entangle- 
ment properties with a thermodynamical potential and 
thus eventually with the partition function. In partic- 
ular one can conclude that, at any temperature (above 
the threshold temperature), no pairwise thermal entan- 
glement exists in the FM (AFM) rings. 

We also determined a simple relation between the mea- 
sure of maximal violation of Bell inequality and the in- 
ternal energy. This in turn allows to explicitly show the 
relation between the concurrence and the violation mea- 
sure. The key ingredient in our derivations has been the 
vast symmetry group of the XXX model. The study of 



the relation between thermal entanglement and thermo- 
dynamical quantities for spin models with a poorer sym- 
metry structure is, we believe, an intriguing challenge for 
future investigations. 
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